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Abstract 

We derive formulas for the matrix elements of the lattice Green function for the discrete Poisson equation 
on an infinite square lattice. The partial difference equation for the matrix elements is solved by reducing 
it to a series of first order difference equations, which can then be solved sequentially. These formulas are 
useful in solving two dimensional Poisson equation problems using the finite difference approximation. 
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I. INTRODUCTION 



In this paper we will derive formulas for the matrix elements of the lattice Green function for 
the discrete Poisson equation on an infinite square lattice. The discrete Poisson equation arises 
when the finite difference approximation is applied to the continuous Poisson equation. The re- 
sults derived here will therefore be of considerable interest for solving any problem that can be 
modeled by a two dimensional Poisson equation [1]. One particularly large area of application is 
in the solution of two dimensional electrostatics and magnetostatics problems. Calculating capac- 
itance, inductance, and charge distributions on conductors with boundary conditions at infinity are 
problems of great practical as well as some theoretical interest fl. 

We will begin by introducing some notation and a formulation of the problem. We let L denote 
the lattice Laplacian operator and then write the discrete Poisson equation as: 

m = \f) (i) 

In the lattice basis we have basis vectors \n) associated with the lattice point r n = n\a\ + n 2 a 2 , 
where a,- = ax.u Xi ■ Xj = S(i, j), and rij is an integer. In this basis eq. [Ubecomes 

Y^UnWn) =f(7l) (2) 
n 

The Laplacian matrix elements are given by 

Li n = -48(r/,r n ) + S(r, + a\ , r n ) + 8(r t - a\ , r n ) + 8(r z + a 2 ,f n ) + 8(F/ - a 2 ,r n ) (3) 
The Green function is defined by the equation LG = —I which in the lattice basis is 

Y^UnGnm = X)£(>7 - r n )G{r n - r m ) = -8(/, m) (4) 

n n 

Using eq. Ogives the following recurrence for the matrix elements of G. 

2 

AG{ri - r m ) - £ [G(r/ + at - r m ) + G(r/ - - f m )\ = 8(Z, m) (5) 

With the notation ?/ — r m = (h — m\)a\ + {h-m^jat = p\a\ + p 2 a 2 , and G(r/ — f m ) = G(p\,p 2 ), 
eq. El becomes 

4G(p l ,p 2 )-G(p l + l,p 2 )-G(p l -l,p 2 )-G(p u p 2 + l)-G(p uP2 -l)=b(p l ,0)8(p 2 ,0) (6) 
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For an infinite lattice, an eigenbasis expansion gives the following formula for the matrix elements 
of gQ 

r( \ 1 F F COSXipiCOSX 2 p2 , , fn , 
G[P1,P2) = / / ~ dxidx 2 (7) 

2K 1 Jo Jo 2 — COSXi — COSJC2 

The problem with this equation is that the integral is divergent for all values of piand p2- We get 
around this problem by using the origin referenced Green function g{p\,pi) = G(0, 0) — G(p\, pi) 
which is then given by the following integral 

1 f* f n 1 -COSJCi/?iCOSX 2j P2 . . , Q . 

g{Pi,P2) = 7r=2 / ~^ dx\dx 2 (8) 

2k z Jo Jo 2 — cos^i — cosx2 

This integral is finite for all finite values of p\ and p 2 - 

g(pi,P2) will also provide a solution to eq. [T]provided that the sum of the source terms, f(r n ) 
over all the lattice sites is equal to zero. To see this, note that the solution to eq. [T]in the lattice 
basis is given by 

<Kr/) = -£G in /(r n ) (9) 

Now if we have 

E/fo)=0 (10) 

n 

Then clearly 

<K?z) =!>«/(?«) (ii) 

II 

Since the coefficients of the recurrence in eq. |5]add to zero we see that g(pi,P2) must obey the 
equation 

4g(pi,P2) ~g{p\ + l,P2) ~g(Pi - 1,P2) ~g(PhP2+ 1) ~g(PhP2 - 1) = -S(pi,0)8(p2,0) 

(12) 

We will now show how to solve this equation for g(pi,p2)- 

II. GENERAL SOLUTION OF THE DIFFERENCE EQUATION 

The key to solving eq. [I2]is the use of the symmetry g{p\iPi) = g(p2,Pi)- This will allow 
eq. [12] to be reduced to a series of first order difference equations for the matrix elements along 



the subdiagonals. The value of the elements along the main diagonal will then provide a unique 
solution to the problem. These values can be directly calculated via eq. [Ell 3]. 
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We begin by letting p\ = p 2 = n in eq. [TJland using the fact that g(p\,pi) = g(p2,Pi), we get 

g(n+l,n) +g(n,n- 1) = 2g(n,n) + ^8(n,0) (13) 

This equation is really just a first order difference equation for the first subdiagonal matrix ele- 
ments. To make this observation explicit we introduce the following notation: ro(n) = g(n, n), and 
ri(n) = g(n, n — 1), r\ (n + 1) = g(n + 1 , n) . Eq. fl3lthen becomes 

n(n + \)+ n (n) = 2r (n) + ^6(n,0) (14) 

This equation for the first subdiagonal matrix elements, r\(n), is easily solved once the diagonal 
matrix elements, ro(n), are known. The initial condition for the equation is n(l) = n(0)- This 
comes from the definition, r\(n) = g(n,n — 1), from which we get ri(l) = g(l,0) and n(0) = 
g(0, — 1 ) = g( 1 , 0) . The value of the diagonal matrix elements was previously determined to be |0] 

= stab (15) 

k—l 



The solution to eq. UHcan of course be found as the sum of the solutions to the two equations: 



n(n+l) + n(n) = 2ro(n) (16) 
ri(n + l)+n(/i) = i8(n,0) (17) 



The solution that comes from eq. [T6jis the complementary solution of eq. [121 The solution 
corresponding to eq. [T7]will have ro(n) = for all n. 

For the next subdiagonal we let p\ = n and p2 = n — 1 in eq. [l2]and we let riin) = g(n, n — 2) 
to get: 

r 2 (n+l) + r 2 (n)=4ri(n)-ro(n)-ro(n-l) (18) 
and in general for the fcth subdiagonal elements we will have 

r k (n + \)+r k (n) = Ar k -\{n) - r k - 2 {n) - r k - 2 (n- 1) (19) 

where r k (n) = g(n,n — k). This equation can easily be solved to give r k {n) in terms of r k -\(n) and 
rk-l{n). 



r,(«) = (-ir i + 1 (r i _ 2 (i-2)+^ 2 ^-l))-^ 2 (»-l)-45;(-l)Vl(»-j) (20) 

7=1 
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We apply this equation iteratively with r\ in) given by either eq. [T6]or eq. [l7]to find a closed form 
solution for r^n). First we look at the case where r\(n) is given by eq. [T7] This leads to the 
following solution: 

r (n) = (21) 
n(n) = (-l) n+1 \ 

r 2 (n) = (-l)"(n-l) 



1 L 2 
4 



At this point we introduce a new function s(n 7 m), defined as 



2/-1 j 




II 


r n-j-l \ 




V j ) 




\k-2j-l I 



S t I n-i\ I n-j + m/2\ 
s(n,m) = £r- 2 J (22) 

In terms of this new function, eq. ETlbecomes 

r k (n) = (-l) n+k ^s(k-l,2(n-k)) (23) 

Now the complementary solution, that starts from eq. [TH can also be expressed in terms of this 
new function as: 

~Y( U n- Ak-l,2(n-j-k) + l) ' S ( k -2j,2(n + j-k)-l] 

> 2/-1 f-T J 2/-1 



r,(n)=r («-/) + (-l) fc - 

71 



i i 2 -/' 1 / i 2 i 

(24) 

where / = |_§J ■ We can verify that these two solutions for A>(n) satisfy eq. |T9]by direct substitution 
and the use of the following identities for the function s(n,m). 

s(n,m) = s(n — l,m) +s(n — l,m+ 1) + s(n,m— 1) (25) 

4s(n,m) = s(n— l,m) +s(n + l,m) — s(n— l,m + 2) — ^(n + l,m — 2) (26) 

Eq. |25]can be derived from the definition of s(n,m) given in eq. |22j Eq. |26]is gotten by applying 
eq. |25]four times. 

Adding the two solutions for r^(n) and using r^(n) = g(n,m), with m = n — k, gives us the 
following equation for the general matrix element of the Green function: 

(-l) m 

g{n,m) = h{l)-\ — s(n — m— 1,2m) 
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(-1) 



m I 



js(n — m— l,2(m — j) + 1) 



+ 



E(-i)' 



2j-l 



n—m 



+ 



(-1) 



m 



I (- 1 ) 



jS(n-m-2j,2(m + j) - 1) 



(27) 



7=1 



27-1 



where / = [n+^±l\ and /*(/) 




III. CONCLUSION 



The symmetry of the square lattice means that we have g(n,m) = g(m,n). This allowed the 
difference equation for g(n, m) to be reduced to a series of first order equations for the subdiagonal 
matrix elements. Solving these equations gave us an expression for the general Green function 
matrix element g(n 7 m). The value of the diagonal matrix elements g{n,n) determine a unique 
solution. This same procedure can also be used to find the lattice Green function for the discrete 
Helmholtz equation as long as the diagonal matrix elements can be calculated directly. An initial 
investigation indicates that it may also be possible to apply this same type of procedure to cubic 
and higher dimensional lattices. For a cubic lattice, the values along the main diagonal g{n,n,n) 
will give a unique solution to the problem, and so on for higher dimensional lattices. Work on 
these problems is currently being pursued. 
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